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We consider the result of acceleration of heavy ions in the slot gap potential of a very 
young pulsar with a hot polar cap. Photodisintegration of the heavy ions in the radiation 



Abstract 



field of the polar cap and pulsar surface gives rise to a flux of energetic neutrons. Some 
fraction of these neutrons interact with target nuclei in the supernova shell to produce 



O 

neutrino and gamma-ray signals which should be observable from very young supernova 
^ ■ remnants in our galaxy for a range of pulsar parameters. 



><J ■ 1 Introduction 



Young supernova remnants are probably sites of acceleration of particles to high energies 
and, as a result of interactions, also sources of high energy neutrinos and 7-rays [0-0. 
Following the occurrence of SN1987A it was expected that high energy 7-rays would be 
detected from this object, but the observations have so far been negative Also, 
some models required the formation of a pulsar during the supernova (SN) explosion 
which has not yet been discovered. Nevertheless, SN1987A renewed interest in predict- 
ing neutrino and 7-ray fluxes from supernova remnants (SNR) for times soon after the 
explosion. For example, the model in which protons are accelerated of at the reverse 
shock which forms in a relativistic wind from a pulsar as a result of its confinement by 
the supernova envelope, proposed by Rees & Gunn fllSf , has been considered in several 
papers [^|, ^|, [R| [171 . In this model, relativistic protons produce 7-rays and neutrinos in 
collisions with the matter in the shell, and observable 7-ray fluxes were predicted for SN 
1987A if the power in relativistic protons was ~ 10 39 erg s . If the pre-supernova star 
had a strong wind, neutrinos and 7-rays could be produced in collisions of relativistic 
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particles accelerated by the supernova shock itself, with the surrounding matter. For this 
scenario, Kirk et al. |1| predicted detectable fluxes of TeV 7-rays for SN 1993J, but not 
for SN 1987A. 

In the present paper, we consider yet another scenario in which particle acceleration 
along open magnetic field lines may take place in the pulsar magnetosphere if a neutron 
star with a high magnetic field and a sufficiently short period forms during the supernova 
explosion. That acceleration of heavy ions by the Crab pulsar may take place, has already 
been suggested, based on modeling the observed synchrotron wisps in the Crab nebula |fT9"[ 



20fl . In these papers it is proposed that a significant fraction of the rotational kinetic 
energy loss of the Crab pulsar is carried by ions which are injected with a rate close to 
the Goldreich & Julian current, and that these ions are accelerated through at least 20% 
of the total voltage available on open field lines. 

It is expected that soon after supernova explosions the surface temperatures of newly 
formed neutron stars are much higher than derived from recent observations of the thermal 



emission from the surfaces of classical pulsars pl|-p3||. Just after formation, a neutron 
star's surface temperature is expected to be very high (~ 10 11 K) but drop very rapidly, 
and models of neutron star cooling predict temperatures of the order of a few 10 6 — 10 
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K during the first ~ 10 years after the explosion [g4|,[25|]. However, the polar cap of a 



neutron star can have a temperature ~ 10 7 K or higher, due to the heating by electrons 



and 7-rays from cascades in the pulsar magnetosphere |26H28| , and their thermal emission 
depends only on the pulsar parameters. 

Nuclei, probably mainly Fe nuclei, extracted from the neutron star surface and ac- 
celerated to high Lorentz factors can be photodisintegrated during propagation through 
the neutron star's radiation field. Relativistic neutrons extracted from Fe nuclei in this 
way propagate straight out and escape from the magnetosphere, interacting in or passing 
through the surrounding supernova remnant shell. Those interacting in the shell produce 
neutrino and 7-ray signals. Those passing through the shell decay into relativistic protons 
which contribute to the pool of Galactic cosmic rays. 

In Sect. 2 we investigate the photodisintegration of nuclei in the radiation field of a 
pulsar. In Sect. 3 we consider the energetics of pulsars and the acceleration of particles 
in the slot gap, as well as polar cap heating. We then compute in Sect. 4 the spectrum 
of neutrons escaping from the pulsar's radiation field, and in Sect. 5 the expected light 
curves and spectra of 7-rays and neutrinos. The acceleration of Fe nuclei in the outer 
gap of classical radio pulsars p 7 !,^, and photodisintegration of Fe on the nonthermal 



radiation of the outer gap is considered in another paper pl |. 



2 Photodisintegration of nuclei in pulsar thermal ra- 
diation 

The possibility of photodisintegration of Fe nuclei during their propagation through ther- 
mal radiation from very young neutron stars has been noted by Bednarek & Karakula 
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The importance of this process can be evaluated by calculating the mean free path, 
Xa(ia, T), for relativistic nuclei with mass number A and Lorentz factor 7^ for photodis- 
integration, typically extracting one nucleon, in black body radiation with temperature 
T which is given by 



where n(e, T) is the photon number density per unit photon energy at photon energy e for 
black body radiation at temperature T, cr^(e*) is the cross section for photodissociation 
of a single nucleon from a nucleus with mass number A |33]| , e* = 7ac(1 — Pa cos 9) is the 
photon energy in the rest frame of the nucleus, [3ac is the velocity of the nucleus, and 
9 is the angle between the direction of the nucleus and a thermal photon. In Fig. [l] we 
show, for illustration, the mean free paths through isotropic black body radiation with 
temperature T, for nuclei with various mass numbers as a function of Lorentz factor. Note 
that the mean free paths computed by us using Eq. |T] differ slightly for low A nuclei in 
comparison with those given by Karakula & Tkaczyk [|33| because the cross section given 
by Eq. 2 of their paper incorrectly separates the contributions from the Breit-Wigner peak 
and the approximately constant cross section at high energies. 

The radiation field we adopt in the present work is anisotropic and the intensity of 
radiation at a particular point in space has the the following properties: in directions 
towards the neutrons star's polar caps the intensity is equal to that of black body radia- 
tion of temperature equal to that of the polar caps; in directions towards other parts of 
the the neutrons star's surface the intensity is equal to that of black body radiation of 
temperature equal to that of the surface; and in all other directions the intensity is zero. 
The radiation field may differ from this because of radiative transfer effects in the neutron 
star's atmosphere due to the presence of a strong magnetic field, and this may cause the 
intensity to depend on the angle to the magnetic pole (see e.g. ref. This effect 

may modify pulse profiles of pulsed X-ray emission, but it would not affect significantly 
the photodisintegration since the modulation with angle varies by a factor less than 4 for 
photon energies ~ 1 keV (typical thermal photon energies for T = 10 7 K). 

To check if photodisintegration of Fe nuclei propagating in the anisotropic radiation 
field of a neutron star perpendicular to its surface is important, we compute the average 
number of neutrons extracted from Fe nuclei as a function of their Lorentz factor 7p e . 
We have made the reasonable approximation, that in every photodisintegration, only one 
nucleon is extracted. In Fig. ^](a) we show the average number of neutrons extracted from 
Fe nuclei as a function of Lorentz factor 7p e for different surface temperatures. It is clear 
that for a neutron star with a uniform surface temperature, the temperature needs to be 
rather high, above ~ 10 7 K, for the process to occur efficiently Such surface temperatures 
are higher than expected except within a few days of the explosion. However, as we shall 
discuss in a later section, neutron star polar caps may have very high temperatures due 
to particle acceleration and cascading in the magnetosphere. We show in Fig. |2|(b) the 
average number of neutrons extracted from Fe nuclei as a function of Lorentz factor in 




(1) 
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the radiation field of a heated polar cap of radius 10 5 cm for three temperatures, and we 
see that efficient photodisintegration can take place in such a radiation field. 

We should point out that the presence of a magnetic field can have an effect on any 
process, including photodisintegration of nuclei by thermal photons, due to a change 
in kinematics, resonant effects, strong field particle wavefunctions, etc. In this paper, 
we negelect such effects which are probably small since heavy nuclei are involved. We 
shall, however, consider in the next section the possibility of photodisintegration of nuclei 
directly on the pulsar's magnetic field. 

3 Photodisintegration of nuclei in the pulsar's mag- 
netic field 

The photodisintegration of iron nuclei may also occur because of the presence of an 
extremely high magnetic field. The component of magnetic field perpendicular to the 
direction of motion of a particle with velocity (3c and Lorentz factor 7, B±(x, 7), at 
distance x = f3ct above the polar cap will give rise to electric and magnetic fields at 
time if = t/7 in the particle's rest frame which are, in SI units, E'±(t') = / ~ycB±(x, , y) 
and B' ± (t') = / ~fB±(x, , j) respectively. These fields are characteristic of electromagnetic 
radiation, and the time-integrated Poynting flux per unit angular frequency is 

rIF' 9 

— \EW)\\ (2) 



daduj' Hqc 

where E'±(uj') is the Fourier transform of E'±(t'), and the number of photons per unit 
area per unit angular frequency is 

" N ' 1 VxMI 2 (3) 



dadoo' hoo 1 (IqC 

This method is referred to as the Weizaeker- Williams method of virtual quanta (see 
Section 15.4 in ref. f35fl). 

As a relativistic iron nucleus travels along a magnetic field line with radius of curvature 
p it will acquire a drift velocity Vd perpendicular to the field direction with order of 
magnitude 

v d ~ — , (4) 

L) B p 

where ujb = 9.6 x 10 7 ZB/A rad s _1 is the non-relativistic cyclotron frequency and B is 
in Tesla (see e.g. ref. p5]| , p 587). Thus, the perpendicular component of magnetic field 
is given approximately by 

B±.{x,i) « — B(x). (5) 
c 
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We shall assume a dipole magnetic field 



B{x) =B s {R/rf (6) 
where r R + x, R is the neutron star's radius, and we take the radius of curvature to 



be that of the last open field line [26 



p = 2 V2 « 1.4 x 10 9 (r 8 P) 1/2 cm (7) 

where r§ = r/10 8 cm, and P is the period in s. 

As a nucleus is extracted from the neutron star's surface and is accelerated along an 
open field line, B± will initially increase as 7 increases and then decrease as p increases 
and B decreases. In the instantaneous rest frame of the nucleus, B' ± (and E'±) will 
also increase from zero to some maximum, and then decrease. In the high u' limit, the 
the spectrum of virtual quanta will be determined by the most rapid change in E'± = 
cB' ± . To see whether the presence of the pulsar's magnetic field is at all important for 
photodisintegration, we consider an extreme case: B± changes instantaneously from zero 
to some constant high value, 

B' ± (X, 7 ) « e(0) 7m ax£ ± (x = 0, 7 = 7max). (8) 

If enhanced photodisintegration due to the magnetic field is unimportant for this case, 
then we can safely neglect it in the present work. 

For the case under consideration, the electric field in the instantaneous rest frame 
varies as 

E' ± (t') « 0(O)P' O (9) 

where 

E' = 7maxCP±(x = 0, 7 = 7 max ). (10) 



Then 



where i = J — 1 and 



dN' 12/ E'. 



o 
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For a photodisintegration threshold of hu' miri ~ 30 MeV (uj' min rs 4.6 x 10 22 rad s ^ we 
obtain an optical depth for photodisintegration on the magnetic field of 



T B ~ CTphoto 



■do;' 



' . daduj 

min 

1 -7/2 ./ 



~ ^photo g^^ ^o^min , (13) 

where <7 p hoto ~ 10~ 26 cm 2 is an approximation to the energy-dependent photodisintegra- 
tion cross section of Fe nuclei above threshold. Taking B s = 1.414 x 10 12 gauss, P = 5 ms, 
R = 1.2 x 10 6 cm, and 7 max = 5 x 10 7 , we obtain E' Q « 5 x 10 19 V m _1 giving r B ps 5 x 1CT 7 . 
Clearly, photodisintegration of Fe nuclei on the pulsar's magnetic field is negligible for the 
presently assumed pulsar parameters. 



4 Pulsar energetics and particle acceleration 

We are interested in the case of a very young pulsar which has just been formed during 
a supernova explosion, and for a few years after formation. It is assumed that pulsar 
acceleration can operate very soon after the formation of a neutron star and can tap into 
a substantial fraction of the spin-down power that is not associated with gravitational 
radiation. We approximate this non-gravitational component of the spin-down power by 
the magnetic dipole radiation power (erg s _1 ) given by 

L cm (tt, B) = B 2 R 6 n 4 sin 2 i/Qc 3 , (14) 

where B is the surface magnetic field at the magnetic pole in G, R is the radius of the 
neutron star in cm, Q = 2tt/P where P is the period in s, and i is the angle between the 
spin axis and rotation axis. 

The total magnetic dipole spin-down power should also include the energy loss due 
to plasma flow (e.g. ref. However, pulsar spin-down through plasma winds is not 



well understood for the case of i ^ 0. Here, we use spin-down through dipole radiation as 
a working hypothesis and B sin i is treated as one parameter (rather than having B and 
sinz separately). 

For the typical magnetic field of a classical pulsar B ~ 10 12 G, for acceleration of ions 
up to sufficiently high energies for photodisintegration to take place we need to have a 
short initial period. Thus, in all the calculations we make in the present paper we shall 
adopt Bsini = 10 12 G, and an initial period of either P = 5 ms or P = 10 ms. 

Equating L cm with the loss of rotational kinetic energy, the period at time t is given 
by the well-known formula, 

P 2 (t) = 1.04 x lO" 15 ^ sin 2 % + P 2 , (15) 

where B = 10 12 i?i2 G, Po is the initial pulsar period in seconds, and we have used 
R = 1.2 x 10 6 cm and I = 1.4 x 10 45 g cm 2 . Note that Eq. [15] neglects other sources 
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of loss of rotational kinetic energy, e.g. by gravitational radiation or electromagnetic 
quadrupole radiation (see the discussion by Manchester and Taylor f38|). The electro- 



magnetic quadrupole radiation which might be important in very young pulsars is very 
uncertain and we shall neglect this. However, gravitational radiation losses are known to 
be important for very young pulsars and the power in gravitational radiation (erg s _1 ) is 
given by 



32 G 



L gv (n, k) = ~—^-i 2 k 2 n b (16) 



where Gat is the gravitational constant, / is the moment of inertia and k is the ellipticity 
p7fl . We adopt k = 3 x 10~ 4 which is required in the case of the Crab pulsar to give 
a pulsar age equal to the time elapsed since the observation of the Crab supernova in 
AD 1054 |37|. 



Starting with an initial pulsar period we have performed numerical integrations over 
pulsar age to obtain the period as a function of age for the case where both electromagnetic 
losses and gravitational radiation losses are included. We show in Fig. |3| the period as a 
function of age for Bsini = 10 12 G and k = 3 x 10~ 4 for two initial periods: Pq = 5 ms 
and Pq = 10 ms. 

According to standard models of neutron star cooling, the surface temperature de- 
creases to T s ~ 10 7 K at a few days after the explosion, and later cools much more slowly, 
reaching T s ~ 4 x 10 6 K in one year. In our calculations, we interpolate/extrapolate the 



results of Nomoto and Tsuruta |24j for which log(T s /l K) is 6.99, 6.72, 6.58, and 6.55 at 
log(t/l y) = -2, -0.94, 0.08, and 1.11 respectively. 

Models of high energy processes in pulsars predict polar cap temperatures can be 
significantly higher because of the heating caused the relativistic e ± pairs and 7-rays 
created in the gaps in the pulsar magnetosphere. For the space charge limited flow 



model P^,H2lH3j a lower limit on the polar cap temperature can be estimated 

Tf « 2.6 x 10 5 p- 19 / 32 Bli 4 Rl /16 K, (17) 

where P is in seconds and R = 10 6 Rq cm. In deriving this formula it is assumed that 
the heating rate (given by Eq. 75 in Arons & Scharlemann |^7f| ) applies to the region of 
the polar cap defined by the last open field line which intersects the pulsar surface at 
(r, 9) = (R, 9 C ) where 

9 C « (fti?/c) 1/2 (18) 

giving a polar cap radius of 

r c ss (2nR 3 /cP) 1/2 . (19) 



Arons & Scharlemann [27] note that the heating rate may actually be even higher, as 
riven by the model of Ruderman & Sutherland [26 , if most of the electrons from pairs 
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created in the polar gap can be reversed and fall on the polar cap region. The Ruderman 
& Sutherland heating rate (given by Eq. 26 in their paper) defines an upper limit to the 
polar cap temperature 

T C RS « 2.8 x 10 6 p- 8 / 28 I^ 2 /28 i? 6 - 17/28 K. (20) 

Note that in both models the polar cap temperatures depends on the pulsar parameters. 
For very rapid pulsars with high magnetic fields they can reach ~ 10 K or higher. In 
the present paper we shall consider three cases: (1) no polar cap heating (the whole star 
has temperature T s ); (2) surface except for polar cap has temperature T s and polar cap 
has the higher of T s or T^ s (Eq. [17]) referred to as "moderate polar cap heating"; and 
(3) surface except for polar cap has temperature T s and polar cap has the higher of T s or 
T* s (Eq. H) referred to as "maximum polar cap heating" . 

For hot polar caps, thermionic emission of ions is allowed, and because of the free 
supply of charges from the polar caps, the rotation-induced potential should be space 
charge limited (e.g. refs. |2"7|, |55|,fl"D"||). Whether ions or electrons are accelerated out along 
field lines depends on the sign of the accelerating electric field. Here we consider only the 
case in which ions are accelerated. 

It is usually assumed that the pulsar surface is composed largely of highly ionized Fe 
although other compositions, e.g. He nuclei |4T| , are also possible. Here, we assume that 
the all the outflowing ions are fully- ionized Fe (e.g. |4(|). Fe nuclei accelerated in the 



magnetosphere may reach sufficiently high Lorentz factors such that photodisintegration 
could occur efficiently. To gain an impression about the pulsar parameters for which this 
may be important, we assume particles are accelerated by the electric field in the slot 



gap f2^, [42], [43|] . In this model the electric potential in the gap is given approximately by 

« 1.13 x 10 2 9 4 C B Rg(x)r](l - rf) sin C sin i V, (21) 

where r\ is the ratio of the colatitude angle of the open field lines at the surface to 9 C , and 
( is the azimuthal angle with respect to the magnetic pole. In the equation, g(x) ~ x 2 for 
x < 0.59 c , and g(x) ~ 9 c [(x + l) a5 — 1] for x > 0.59 c , and x is the distance along the gap 
in units of R. In writing Eq. [H], we ignore the terms due to the complicated geometry 
of the gap. The ion acceleration zone corresponds to i > 7r/2. We shall adopt 77 = 0.1 
and £ = 37r/2 in our work. This potential is shown in Fig. |] for B sini = 10 12 G and for 
periods equal to the two initial periods, i.e. P = 5 ms and P = 10 ms. 



5 Spectrum of neutrons from photodisintegration of 
Fe nuclei 

The extraction of energetic neutrons from accelerated nuclei has important consequences 
for the transport of energy from the pulsar's vicinity to the nebula since they move bal- 
istically through the magnetosphere and beyond the light cylinder. Their Lorentz factors 
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are approximately equal to the Lorentz factors of the parent nuclei, and are sufficiently 
high that they typically decay at large distances from the pulsar. Initially the SNR shell 
is opaque to neutrons so that most neutrons interact in the dense shell. Later, when the 
shell becomes more transparent to neutrons, they typically move through the SNR shell 
and decay outside it. Most of the energy of the neutrons is then taken by protons which 
are trapped and isotropized by the ambient magnetic field, and these protons then wait 
for target nuclei, i.e. the SNR shell, to arrive. When the SNR shell catches up with the 
relativistic protons, production of 7-rays and neutrinos occurs through hadronic collisions 
but, as we shall discuss in Section 6, this contribution is negligible for the present radia- 
tion field and electric potential. The protons will also contribute to the Galactic cosmic 
rays, although their contribution is relatively small. 

The dependence of the electric potential in the case of the slot gap model as a function 
of distance (Eq. |2l| ) defines the Lorentz factor of Fe nuclei during their propagation from 
the surface. Using this potential we compute the average number of neutrons extracted 
from Fe nuclei as a function of pulsar age. We show this in Fig. |5|(a) for P = 5 ms for each 
of the three radiation field models: no polar cap heating, moderate polar cap heating, 
and maximum polar cap heating. Note that with maximum polar cap heating the polar 
cap temperature is not significantly higher than the stellar surface for t < 2.2 days, and 
with moderate polar cap heating the polar cap temperature is not significantly higher 
than the stellar surface for t < 18 days. In Fig. H(b) we show the average cumulative 
number of neutrons extracted from a single Fe-nucleus during acceleration for P = 5 ms 
as a function of distance along the slot gap for the three radiation field models. 

Since at a particular distance, a neutron extracted from the accelerated ion will have 
an energy determined by the slot gap potential, we can now compute the energy spectra 
of neutrons escaping from the pulsar radiation field per single Fe nucleus accelerated. 
In Fig. ^| we show the spectrum of neutrons extracted from a single Fe-nucleus during 
acceleration, N n (E n ), for a pulsar with Pq = 5 ms and Bsini = 10 12 G, for the three 
radiation field models. Results are given for (a) t = 1.15 days (P = 5 ms), (b) t = 
20.5 days (P = 5.05 ms), and (c) t — 1 year (P = 5.69 ms). In Fig. [7| we compare the 
spectrum of neutrons at t — 1 year produced for P = 5 ms with that produced with 
P = 10 ms. 

The differential production rate, N n (E n ), of neutrons extracted from Fe nuclei depends 
on the total rate of Fe nuclei injected, N Fe , and is given by 



where L em (Q, B) is the magnetic dipole radiation approximation to the total non-gravitational 
radiation spin-down power (Eq. [14]), $ max (f2, B) = $(r LC /i?) is the maximum acceleration 
potential traversed by Fe nuclei given by Eq. ^l| and tlc is the distance to the light cylin- 
der up to which this acceleration is assumed to occur, Z = 26 is the atomic number of Fe, 
and e unit electric charge. The parameter £ gives the fraction of total pulsar power used 
to accelerate Fe nuclei and has an upper-limit £ « $ max /$o ~ 0.29 where $0 = 0.56^BR 



N n (E n ) = N Fe N n (E n ) 




N n (E n ) 



(22) 



Ze$ max (n,£?) 
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is the maximum potential drop across the polar cap (e.g. Luo et al. Note that 

we write the rate of injection of iron nuclei as the ratio of a fraction of the power, £L em 
to the maximum particle energy, Ze$ max since the ratio gives the same order-magnitude 
estimate of the injection rate as that derived from Goldreich- Julian density. 



6 Gamma-ray and neutrino spectra 

We make the approximation that in shell-type SNR the shell has a thickness which is much 
smaller than the radius of the shell, r$N- The optical depth of the shell to nucleon-nucleon 
collisions is then 



T, 



PP 



M SN a pp /(A7rr 2 SN m p ) « 3.5 x l^M w /{f3 m t y Y (23) 



where M$n = lOMio M Q is the mass ejected into the shell during the explosion, a. 
3 x 10 -26 cm 2 is the proton-proton inelastic cross section, m p is the proton mass, t = t y 
years is the time after the explosion, and t>sN = Psnc is the shell expansion velocity. 
Initially the shell is optically thick. For example, if M$n — 1OM and /3sn — 0.03 the 
optical depth is greater than 1 for t < 2 months. All the results we shall present below 
are for /? SN = 0.03 and M SN = 10 M . 

During the early phase most of the neutrons interact with matter in the shell producing 
7-ray and neutrino fluxes. These signals may be observed provided the 7-ray and neutrino 
beams intersect the direction to the Earth, or sweep across the Earth giving rise to a pulsed 
signal. For a beaming solid angle of fib steradians, the neutrino flux may be calculated 
from 

F V {E V ) « ^[1 - exp(-r pp )] J N n (E n )P^(E u , E n ) dE n (24) 

where d is the distance to the SNR, and P£,(E U , E n )dE u is the number of neutrinos 
produced with energies in the range E u to (E u + dE u ) (via pion production and subsequent 
decays) in multiple nucleon-nucleon interactions of a nucleon of energy E n . Nucleon- 



nucleon interactions were treated as described in |45[ and the pions produced were decayed 
using SIBYLL [^BJ. All the results given below are for fij, = 1 sr, and for smaller beams 
the fluxes should therefore be scaled up by a factor . 

In the case of 7-ray production, a significant fraction of the 7-rays will interact by 
pair production with target nuclei in the shell and will not be observed. Also, during the 
first few months, TeV 7-rays could interact with radiation from the photosphere of the 
supernova [}i7|| . By several days after the explosion, the radius of the photosphere becomes 



smaller than that of the SNR shell, where the neutrons interact, and continues to decrease 
in radius relative to the shell. Thus, interactions with radiation from the photosphere is 
likely to be a small effect here as the 7-rays will be traveling in the directions of their 
parent neutrons, i.e. radially outwards from the neutron star, and the angle between the 
directions of the 7-rays and soft photons from the photosphere will be small, reducing 
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the interaction probability. Concentrating on pair production with matter, the optical 
depth of the shell is determined by the mean free path for pair production which at high 
energies is (9/7) of the radiation length, giving r 7P « 0.7t pp . For a neutron interacting 
at some fraction, /, of the way through the shell, a 7-ray produced at that point would 
have a probability of exp[— (1 — /)r 7P ] of escaping from the shell. Integrating over neutron 
interaction points (7-ray emission points) within the shell we arrive at the gamma ray 
flux at Earth, 

W ^ 6XP( ^j r ~) X _ P( r T7p) / N n (E n )P%(E 7 , E n ) dE n , (25) 

where P^£(Ej, E n )dE 7 is the number of 7-rays produced with energies in the range E™ 
to {E-y + d£/ 7 ) (via pion production and subsequent decays) by multiple nucleon-nucleon 
interactions of a nucleon of energy E n . Note that this neglects cascading in the matter as 
a result of bremsstrahlung and subsequent pair production, and so the 7-ray flux at low 
energies will be somewhat underestimated. 

Light curves for neutrinos above 1 TeV, and 7-rays above 100 MeV and 1 TeV, are 
shown for Bsini = 10 12 G in Fig. |j|(a) for P Q = 5 ms and in Fig. |j(b) for P = 10 ms. 
The results are shown here, and in the remaining figures, for a distance of d = 10 kpc. 
In both cases, results are given for three radiation field models (no polar cap heating, 
moderate polar cap heating, and maximum polar cap heating), and the pulsar period and 
the surface and polar cap temperatures vary appropriately with time after the explosion. 
In Fig. ||(a) we show the energy spectrum of 7-rays, and in Fig. [|(b) the energy spectrum 
of neutrinos, both at t = 0.1 year for each initial period and radiation field model. In 



Fig. [lCj we show how the energy spectrum of neutrinos evolves with time for the case of 

we 



Po = 5 ms and the maximum polar cap heating model. In both Figs. [|(b) and 10 
show the atmospheric neutrino background flux within 1° and within 10° of the source 
direction based on the intensity calculated by Lipari |48| . 



7 Discussion and Conclusion 

As we shall see, we require fairly high acceleration efficiencies in very young pulsars for 
observable fluxes of 7-rays and neutrinos. However, high energy 7-ray observations of 
several pulsars appear to suggest an increase in 7-ray efficiency with age. This may 
be attributed to the relation between acceleration efficiency for primary particles that 
start the electron-positron pair cascade and the age, provided the pulsed 7-ray emission 
from pulsars is from pair cascades. This may be true for young pulsars with ages more 
than ~ 10 3 years. For very young pulsars (a few months to a few years) as discussed in 
this paper, thermal emission from the polar caps or neutron star's surface may have an 
important effect on particle acceleration. For very young pulsars with strong magnetic 
fields and hot polar caps, the efficiency may be higher for younger pulsars. This is because 
energy losses due to resonant inverse Compton scattering are important and can prevent 
primary electrons or positrons from starting a cascade in regions near the polar caps [[19[ 
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We next discuss whether the signals we predict from young SNR are observable with 
existing and planned experiments (see e.g. ref. [|5(J for a summary of 7-ray detector 
thresholds). The 7-ray light curve peaks at about 2 months after the explosion (see 
Fig. ||). The peak flux above 100 MeV for a pulsar at 10 kpc with Psinz = 10 12 G 
ranges from 2 x 10~ 10 cm~ 2 s _1 (Po = 10 ms, no polar cap heating) to 4 x 10~ 9 cm -2 s _1 
(Po — 5 ms, maximum polar cap heating) if ^fi^ 1 = 1 sr -1 . The sensitivity of the EGRET 
detector on the Compton Gamma Ray Observatory for 100 MeV 7-rays is ~ 7 x 10~ 8 
cm -2 s _1 , and with this sensitivity 7-rays from a Galactic source should easily be detected 
for several months if a significant fraction of the electromagnetic power (£ ~ 1) goes into 
accelerating heavy nuclei and the pulsar's 7-ray beam sweeps across the Earth and has 
Qb < 0.07 sr (Po — 5 ms, maximum polar cap heating) or Qb < 3 x 10~ 3 sr (Po = 10 ms, no 
polar cap heating). With future 100 MeV 7-ray detectors, even lower detection thresholds 
(e.g. ~ 5 x 10 -9 cm~ 2 s _1 for GLAST) should enable detection even if Qb were larger. 

The current generation of Cherenkov telescopes operating at TeV energies (e.g. the 
Whipple Observatory, and the CANGAROO telescope [|51jj) have thresholds of a ~ 3 x 
10~ 12 cm -2 s , and GRANITE III (Whipple telescope with new camera) is expected to 
have a threshold of ~ 10 -12 cm~ 2 s _1 . For Galactic sources, the predicted TeV fluxes (see 
Fig. H) are well above the Whipple threshold and should easily be detected even for the 
least optimistic case (Po = 10 ms, no polar cap heating) provided < 3.8 sr, for the 

most optimistic case (Po = 5 ms, maximum polar cap heating) provided £,~ l £lb < 75 sr. 
With the sensitivity of current TeV 7-ray telescopes, it should be possible to detect 7- 
rays from sources in the Magellanic Clouds provided £,~ l flb is sufficiently high, and we 
shall discuss below constraints placed on our model by the non-detection of 7-rays from 
SN 1987A. 

Available upper limits on the TeV 7-ray flux from SN 1987A are as follows: t ~ 10 
months (November 1987) — 2.3 x 10~ n cm -2 s" 1 above 1 TeV (Raubenheimer et al. ||12|| ); 
t 1 year — 6.1 x 10~ 12 cm" 2 above 3 TeV (Bond et al. |§); t ~ 1 year (January/February 
1988) — 1.6 x 10~ 10 cm" 2 s _1 above 0.4 TeV (Chadwick et al. JTJ). Assuming a distance of 
55 kpc, and if the observer is inside the beam of energetic neutrons injected by the pulsar, 
which will also define the 7-ray beam, then the value of can be constrained by the 

1 TeV upper limit mentioned above, giving > 0.36 sr for the case of Po = 5 ms 

(see Fig. |8]a). However, we note that since no pulsar has yet been discovered inside 
SN 1987A this could be interpreted as suggesting that we are outside the beam of the 
neutron injection or that no pulsar was formed. 

The technique for constructing a large area (in excess of 10 4 m 2 ) neutrino telescope 



has been known for more than a decade [ 52| . In pioneering work, the DUMAND Col 



laboration developed techniques to instrument a large volume of water in a deep ocean 
trench with strings of photomultipliers to detect Cherenkov light from neutrino-induced 
muons. Locations deep in the ocean shield the detectors from cosmic ray muons. The 
DUMAND detector |kJ was designed to be most sensitive to neutrinos above about 1 
TeV, and prototypes of other similar experiments are already in operation such as that in 



Lake Baikal, Siberia [54], and NESTOR off the coast of Greece [55]. An exciting recent 
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development has been the construction of a DUMAND type detector deep in the polar 
ice cap at the South Pole. This experiment called AMANDA uses the same principle as 
DUMAND but takes advantage of excellent transparency of the polar ice under extreme 
pressures and a stable environment in which to embed the detectors f56|. These exper- 



iments which operate typically above 1 TeV have the potential to be expanded in the 
future to a detector on the 1 km 3 scale. The background to these experiments at 1 TeV 
is due primarily to atmospheric neutrinos, but at higher energies there is an uncertain 



background of prompt muons from charm production (see e.g. Gaisser et al. |57[ for a 
discussion) . For recent reviews of the predicted intensity of diffuse high energy neutrinos 
of astrophysical origin see refs. |]57|,p8| 



Referring to Figs. |9|(b) and [K] we see that for good angular resolution (1°) the neutrino 
flux above 1 TeV is well above the atmospheric background, and even for 10° resolution 
the neutrino signal may be observable. However, observation of any signal would require 
the neutrino beam (essentially the beam of the primary neutrons) to sweep across the 
Earth. For 1° resolution and ^fi^ 1 = 1 sr _1 , TeV neutrinos should be easily observable if 
their flux remained constant at the level given in Fig. |S|(b) for 0.1 y. However, the neutrino 
light curve drops rapidly at t > 0.1 y (see Fig. |8|(b)) and so very large neutrino telescopes 
will be required. With 10° resolution, detection of galactic sources at 1 TeV would be 
marginal even in the most optimistic scenario except for the case of extremely large 
detectors, e.g. km 3 [|59| which, with large statistics, could see a weak signal significantly 
below the atmospheric background. Note, however, that the emission peaks at energies 
~ 10-100 TeV, and that telescopes sensitive in this range should be able to detect a 
neutrino signal for the P = 5 ms case from a galactic source even if the resolution were 
10°. A large detector with excellent angular resolution would be required to observe very 
young SNR in the Magellanic Clouds unless the emission were beamed strongly towards 
the Earth. 

We have neglected any signals arising from the protons extracted from Fe-nuclei, and 
also signals from nuclei that were not completely fragmented into nucleons in the radiation 
field of the pulsar. We assume these particles, being charged, will be trapped in the central 
region of the SNR where the matter density is expected to be relatively low, and would 
not contribute significantly to the 7-ray and neutrino fluxes. These nuclei would be 
subject to adiabatic deceleration as the SNR expands. However, some would accumulate 
and eventually contribute to the galactic cosmic rays towards the end of the SNR's life. 
Neutrons that do not interact promptly travel far from the pulsar where they decay into 
protons which await the arrival of target nuclei in the supernova shell, and then produce 
a delayed isotropic neutrino and 7-ray signal. We have estimated these delayed neutrino 
and 7-ray fluxes for the present radiation field models and slot gap potential. The 7-ray 
light curves peak at about 6 months with fluxes about 10~ 3 of the peak flux from neutron 
collisions in the shell for i\ — 1 sr. Thus, unlike the case of 7-rays from interactions of 
neutrons with the shell, the delayed fluxes are unlikely to be observable. The TeV neutrino 
light curve is of the same magnitude as the TeV 7-ray light curve. Those protons from 
neutron decay outside the nebula also contribute to the Galactic cosmic rays. We estimate 
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the total energy going into cosmic rays, with typical energies of 10 6 GeV, integrated over 
the pulsar's age to be ~ 10 45 erg, which is much less than the ~ 10 50 erg potentially 
available through shock acceleration. 

In the case of filled nebulae, the 7-ray and neutrino fluxes could be higher than esti- 
mated here. For example, t pp would be a factor of 3 higher than given by Eq. ^ in the 
case of uniform filling. More importantly, cosmic rays trapped inside the nebula would 
also continuously contribute to the 7-ray and neutrino fluxes. Such a model is discussed 



in the context of the Crab Nebula in another paper [31 



We have investigated the use of higher magnetic fields, e.g 10 13 G, and found that this 
has two main effects: (1) the polar cap heating is increased resulting in more neutrons 
extracted per Fe nucleus (but up to a maximum of 30), (2) the maximum acceleration 
potential is increased. The second effect causes the peak of the 7-ray and neutrino emission 



to move to higher energies, and the injection rate of Fe nuclei to decrease (see Eq. 22 



Overall, the net effect of using extremely high magnetic fields is unexpectedly to reduce 
the predicted fluxes to levels comparable to or lower than those predicted here for typical 
pulsar magnetic fields. 

We should point out that in this paper we have ignored the fall-back of matter onto 
the neutron star. If fall-back occurs, which cannot be ruled out, there are two immediate 
consequences: acceleration of lighter ions such as carbon nuclei may occur, or acceleration 
is completely quenched. In the present work we have also ignored production of 7-ray 
lines. The neutrons and protons from the Fe photodisintegration will excite nuclei of 
C, N, O, etc. in the supernova shell to make gamma-ray lines in the 1-10 MeV range. 
These lines may be detectable after about 1 year, when the shell becomes optically thin 
to scattering and absorption. Calculation of the flux of these lines is however beyond the 
scope of the present paper. 

In Conclusion, we have considered a model for 7-ray and neutrino emission by very 
young supernova remnants in which ions are accelerated in the slot gap of a highly magne- 
tized rapidly spinning pulsar. Energetic neutrons are extracted from the ions by photodis- 
integration during interactions with thermal radiation from the neutron star surface and 
hot polar caps. Gamma ray and neutrino signals are produced by energetic neutrons in- 
teracting with target nuclei as they travel out through the SNR shell. For a limited range 
of pulsar parameters these signals should be observable from very young SNR in our 
Galaxy with existing and planned 7-ray and neutrino telescopes within one or two years 
of the supernova explosion. The emission peaks above TeV energies and so ground-based 
optical Cherenkov detectors are more sensitive than the lower energy satellite telescopes. 
Neutrino telescopes should also be able to detect the predicted neutrino signals which 
peak at ~ 10 TeV energies. 
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2 3 4 5 6 

Figure 1: Mean free path for photodisintegration of nuclei in isotropic black body radiation 
at temperature T = 10 7 K as a function of Lorentz factor, 7^4, for various mass numbers 
A: 4 (top curve), 8, 16, 40, and 56 (bottom curve). 
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Figure 2: (a) Average number of neutrons extracted from a single Fe-nucleus during 
propagation along the neutron star's magnetic axis, from its surface at the pole to infinity, 
as a function of Lorentz factor. Neutron star spherical surface of radius 1.2 x 10 6 cm at 
constant temperature T = 5 x 10 6 K (bottom curve), 10 7 K, 1.5 x 10 7 K, and 2 x 
10 7 K (top curve), (b) Average number of neutrons extracted from a single Fe-nucleus 
during propagation along the neutron star's magnetic axis, from its surface at the pole 
to infinity, as a function of Lorentz factor, for field of heated polar cap of radius 10 5 cm 
with temperature T = 10 7 K (bottom curve), 2 x 10 7 K, and 3 x 10 7 K (top curve). 
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Figure 3: Evolution of pulsar period for Bsini = 10 12 G for P = 5 ms (lower curves) 
and P — 10 ms (upper curves). Dotted curves are for electromagnetic radiation losses 
only, and solid curves include the effect of gravitational radiation losses. 
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Figure 4: Slot gap potential given by Eq. ^T] for Bsini = 10 12 G and P = 10 ms (lower 
curve) and P = 5 ms (upper curve). The potential is assumed constant for x larger than 
the light cylinder radius. 
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Figure 5: Average cumulative number of neutrons extracted from a single Fe- nucleus 
during acceleration by the assumed slot gap potential for Bsini = 10 12 G and P = 5 
ms through the pulsar's radiation field: (a) as a function of neutron star age; (b) as a 
function of distance along the slot gap at t = 20.5 days. Results are shown for no polar 
cap heating (solid curve), moderate polar cap heating (dotted curve), and maximum polar 
cap heating (dashed curve). 
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Figure 6: The spectrum of neutrons extracted from a single Fe-nucleus during acceleration 
for a pulsar with Pq = 5 ms and Bsini = 10 12 G at (a) t — 1.15 days (P = 5.00 ms), (b) 
t = 20.5 days (P = 5.05 ms), and (c) t — 1 year (P = 5.69 ms). Results are shown for the 
three radiation field models: no polar cap heating - solid histogram; moderate heating - 
dotted histogram; maximum heating - dashed histogram. 
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Figure 7: The spectrum of neutrons extracted from a single Fe-nucleus during acceleration 
at t — 1 year and Bsini = 10 12 G for the following initial pulsar periods Pq = 5 ms 
(thin histograms), and Pq = 10 ms (thick histograms). Results are shown for two of 
the radiation field models: no polar cap heating - solid histograms; maximum heating - 
dashed histograms. 
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Figure 8: The dependence of the integral flux of neutrinos [y^ + v^) above 1 TeV (solid, 
dot-dash, and dot-dot-dot-dash curves), 7-rays above 100 MeV (dotted curves), and 7- 
rays above 1 TeV (dashed curves) as a function of time t measured from the explosion for 
d — 10 kpc and _Bsini = 10 12 G and for (a) P = 5 ms, and (b) P = 10 ms. Results are 
shown for the three radiation field models: no polar cap heating - dot-dot-dot-dash curve 
and lowest dotted and dashed curves; moderate heating - dot-dash curve and middle 
dotted and dashed curves; maximum heating - solid curve and upper dotted and dashed 
curves. 
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Figure 9: The spectra for d — 10 kpc at t = 0.1 year of (a) 7-rays and (b) neutrinos 
(Ufj, + u^) produced by collisions of neutrons with the matter of the supernova shell. In 
each case, the upper three histograms are for Pq = 5 ms, and the lower three histograms 
for Pq = 10 ms. Results are shown for B sin i = 10 12 G for the three radiation field models: 
no polar cap heating - lower histogram; moderate heating - middle histogram; maximum 
heating - upper histogram. The atmospheric neutrino background flux multiplied by E 2 
(GeV cm~ 2 s _1 ) within 1° and within 10° of the source direction, based on the intensity 
calculated by Lipari is shown by the hatched bands (each band shows the range of 
atmosheric neutrino background as the zenith angle changes). 
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Figure 10: Evolution of the spectrum of neutrinos (y^ + z/ M ) produced by collisions of 
neutrons with matter in the supernova shell for Bsmi = 10 12 G, Pq = 5 ms and the 
maximum polar cap heating model. Spectra are shown for d = 10 kpc and for t = 10~ 4 y 
(solid histogram); 10~ 3 y (dotted histo gram); 10 2 y (dashed histogram); 10 1 y (dot- 
dash histogram); 1 y (dot-dot-dot-dash histogram). See Fig. ^ for explanation of hatched 
bands. 
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